home *** CD-ROM | disk | FTP | other *** search
/ 8bitfiles.net/archives / archives.tar / archives / commodore-users-of-norman / CUON_121_(06-1985).d64 / planets ra_dec (.txt) < prev    next >
Commodore BASIC  |  2019-04-13  |  4KB  |  159 lines
  1. 10 REM -- ADAPTED FROM 'CELESTIAL BASIC', PG. 78, BY ERIC BURGESS
  2. 11 :
  3. 20 PRINT"[147]"TAB(17)"RADEC"
  4. 30 PRINT"   R.A. & DECLINATION FOR THE PLANETS":FOR I=1 TO 3000:NEXT
  5. 35 :
  6. 36 REM ----------------
  7. 37 REM -- INITIALIZE --
  8. 38 REM ----------------
  9. 39 :
  10. 40 DEF FN ASN(X) = ATN(X/SQR(-X*X+1))
  11. 50 DEF FN ACO(X) = -ATN(X/SQR(-X*X+1))+1.5707963
  12. 60 DEF FN RAD(X) = .01745328 * X
  13. 70 DEF FN DEG(X) = 57.29578 * X
  14. 80 DIM P$(8), PD(8,8), A(8), D(8), L(8), Q(8), R(8), V(8)
  15. 90 K1=6.28318 :REM  CONSTANT USED IN SUBROUTINE FOR HELIOCENTRIC LONG (A)
  16. 175 :
  17. 176 REM ----------------------
  18. 177 REM -- READ PLANET DATA --
  19. 178 REM ----------------------
  20. 179 :
  21. 190 FOR I=0 TO 8: READ P$(I)
  22. 200 FOR J=0 TO 8: READ PD(I,J)
  23. 210 NEXT J,I
  24. 217 :
  25. 218 REM -- PLANET DATA: EPOCH 1960,1,1
  26. 219 :
  27. 220 DATA MERCURY
  28. 222 DATA .071422,3.8484,.388301,1.34041,.3871,.07974,2.73514,.122173,.836013
  29. 224 DATA VENUS
  30. 226 DATA .027962,3.02812,.013195,2.28638,.7233,.00506,3.85017,.059341,1.33168
  31. 228 DATA EARTH
  32. 230 DATA .017202,1.74022,.032044,1.78547,1,.017,3.33939,0,0
  33. 232 DATA MARS
  34. 234 DATA .009146,4.51234,.175301,5.85209,1.5237,.141704,1.04656,.03142,.858702
  35. 236 DATA JUPITER
  36. 238 DATA .00145,4.53364,.090478,.23911,5.2028,.249374,1.76188,.01972,1.74533
  37. 240 DATA "SATURN"
  38. 242 DATA .000584,4.89884,.105558,1.61094,9.5385,.534156,3.1257,.043633,1.977458
  39. 244 DATA URANUS
  40. 246 DATA .000205,2.46615,.088593,2.96706,19.182,.901554,4.49084,.01396,1.28805
  41. 248 DATA NEPTUNE
  42. 250 DATA .000104,3.78556,.016965,.773181,30.06,.27054,2.33498,.031416,2.29162
  43. 252 DATA PLUTO
  44. 254 DATA .000069,3.16948,.471239,3.91303,39.44,9.86,5.23114,.300197,1.91812
  45. 280 :
  46. 281 REM ------------------
  47. 282 REM -- MAIN PROGRAM --
  48. 283 REM ------------------
  49. 287 :
  50. 288 REM -- ENTER DATE
  51. 289 :
  52. 290 PRINT"[147]ENTER THE DATE:"
  53. 300 INPUT"YEAR";Y
  54. 310 IF Y<50  THEN Y=Y+2000:PRINTY"ASSUMED"
  55. 320 IF Y=>50 AND Y<100 THEN Y=Y+1900: PRINTY"ASSUMED"
  56. 330 IF Y<1800 OR Y>2100 THEN GOSUB 50000:GOTO 300
  57. 370 INPUT"MONTH (1-12)";M:IF M<1 OR M>12 THEN GOSUB 50010:GOTO 370
  58. 380 INPUT"DAY";D
  59. 390 IF D>31 AND (M=1ORM=3ORM=5ORM=7ORM=8ORM=10ORM=12) THENGOSUB50010:GOTO380
  60. 400 IF D>30 AND (M=4 OR M=6 OR M=9 OR M=11) THEN GOSUB 50010:GOTO 380
  61. 410 IF(D>29ANDM=2)OR(D=29ANDM=2ANDY/4<>INT(Y/4))THEN GOSUB 50010:GOTO 380
  62. 419 :
  63. 420 REM -- CALC. GREGORIAN DAYS FROM
  64. 430 REM -- EPOCH 1960,1,1 TO DATE
  65. 431 :
  66. 440 DG=365*Y+D+((M-1)*31)
  67. 450 IF M>=3 THEN 490
  68. 460 REM -- CALC FOR JAN & FEB
  69. 470 DG=DG+INT((Y-1)/4)-INT((.75)*INT((Y-1)/100+1))
  70. 480 GOTO 510
  71. 490 REM -- CALC FOR MAR THRU DEC
  72. 500 DG=DG-INT(M*.4+2.3)+INT(Y/4)-INT((.75)*INT((Y/100)+1))
  73. 510 NI=DG-715875
  74. 799 :
  75. 800 REM -- CALC, PRINT INFO FOR PLANETS
  76. 801 :
  77. 820 PRINT"[147]PLANETARY DATA FOR Y/M/D:"Y;M;D
  78. 830 PRINT"WHICH IS ";NI;"DAYS FROM EPOCH 1960
  79. 840 [153]"LENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLEN"
  80. 850 [153]"PLANET  HELIO     DIST     R.A.    DEC"
  81. 860 [153]"         LONG   TO PLANET  HRS     DEG"
  82. 870 [153]"--------------------------------------"
  83. 880 [153]
  84. 890 [129] J[178]0 [164] 8:[141] 1140
  85. 900 A(J)[178]A:D(J)[178]D:L(J)[178]L
  86. 915 [130]
  87. 920 [129] I[178]0 [164] 8
  88. 930 [143] -- SKIP EARTH
  89. 940 [139] I[178]2 [167] [130]
  90. 950 [141] 1280
  91. 960 Q(I)[178]Q:R(I)[178]R:V(I)[178]V
  92. 970 [130]
  93. 980 [129] I[178]0 [164] 8:A(I)[178][165] DEG(A(I))
  94. 985 [143] -- SKIP EARTH
  95. 990 [139] I[178]2 [167] [130]
  96. 995 [153] P$(I);
  97. 1000 [153] [163]8)[197]([202]([196](A(I)),2,5));
  98. 1010 [153] [163]16)[197]([202]([196](Q(I)),2,5));
  99. 1020 [153] [163]25)[197]([202]([196](R(I)),2,5));
  100. 1030 [153] [163]32)[197]([200]([196](V(I)),5))
  101. 1035 [153]
  102. 1040 [130]
  103. 1045 [153]"ONLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLENLEN"
  104. 1050 [153]"AGAIN? (Y/N)":[151]198,0
  105. 1060 [161]A$:[139]A$[178]""[167]1060
  106. 1070 [139] A$[178]"Y" [167] 290
  107. 1080 [139] A$[178]"N" [167] [153]"LOAD":[128]
  108. 1090 [137] 1050
  109. 1100 [137]1100
  110. 1118 :
  111. 1119 [143] -----------------------------
  112. 1120 [143] -- CALCULATION SUBROUTINES --
  113. 1121 [143] -----------------------------
  114. 1122 :
  115. 1125 [143] -- CALC A,D, AND L:
  116. 1126 :
  117. 1130 [143] -- HELIOCENTRIC LONG, A
  118. 1140 A[178]NI[172]PD(J,0)[170]PD(J,1)
  119. 1150 [139] A[177]K1 [167] A[178]((A[173]K1)[171][181](A[173]K1))[172]K1
  120. 1160 [139] A[179]0 [167] A[178]A[170]K1: [137] 1160
  121. 1170 C[178]PD(J,2)[172][191](A[171]PD(J,3))
  122. 1180 A[178]A[170]C
  123. 1190 [139] A[177]K1 [167] A[178]A[171]K1
  124. 1200 [139] A[179]0 [167] A[178]A[170]K1: [137] 1200
  125. 1210 [143] -- DIST PLANET TO SUN, D
  126. 1220 D[178]PD(J,4)[170]PD(J,5)[172][191](A[171]PD(J,6))
  127. 1230 [143] -- DIST PLANET TO ECLIPTIC, L
  128. 1240 L[178]PD(J,7)[172][191](A[171]PD(J,8))
  129. 1250 [142]
  130. 1255 :
  131. 1260 [143] -- CALC R, V, AND Q
  132. 1270 :
  133. 1280 Z[178]A(2)[171]A(I)
  134. 1290 [139] [182](Z)[177][255] [175] Z[179]0 [167] Z[178]Z[170]K1
  135. 1300 [139] [182](Z)[177][255] [175] Z[177]0 [167] Z[178]Z[171]K1
  136. 1310 [143] -- DIST PLANET TO EARTH, Q
  137. 1320 Q[178][186](D(I)[174]2[170]D(2)[174]2[171]2[172]D(I)[172]D(2)[172][190](Z))
  138. 1330 [143] -- ANGULAR DIST FROM SUN, X
  139. 1340 P[178](D(I)[170]D(2)[170]Q)[173]2
  140. 1350 X[178]2 [172] [165] ACO([186](((P[172](P[171]D(I)))[173](D(2)[172]Q))))
  141. 1360 [143] -- RIGHT ASCENSION, R
  142. 1370 [139] Z[179]0 [167] R[178][165] DEG(A(2)[170] [255][171]X)[173]15
  143. 1380 [139] Z[177]0 [167] R[178][165] DEG(A(2)[170] [255][170]X)[173]15
  144. 1390 [139] R[177]24  [167] R[178]R[171]24: [137] 1390
  145. 1410 [139] R[179]0   [167] R[178]R[170]24: [137] 1410
  146. 1420 [143] -- DECLINATION, V
  147. 1430 [139] Z[179]0 [167] V[178][191](A(2)[170] [255][171]X)[172]23.44194[170] [165] DEG(L(I))
  148. 1440 [139] Z[177]0 [167] V[178][191](A(2)[170] [255][170]X)[172]23.44194[170] [165] DEG(L(I))
  149. 1450 X[178][165] DEG(X)
  150. 1460 [142]
  151. 40000 [144]
  152. 49995 :
  153. 49996 [143] --------------------
  154. 49997 [143] -- ERROR MESSAGES --
  155. 49998 [143] --------------------
  156. 49999 :
  157. 50000 [153]"PLEASE USE A YEAR BETWEEN 1800 AND 2100":[142]
  158. 50010 [153]"OUT OF RANGE":[142]
  159.